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We study the free decay of capillary turbulence on the charged surface of liquid hydrogen. We find 
that decay begins from the high frequency end of the spectral range, while most of the energy remains 
localized at low frequencies. The apparent discrepancy with the self-similar theory of nonstationary 
wave turbulent processes is accounted for in terms of a quasi-adiabatic decay wherein fast nonlinear 
wave interactions redistribute energy between frequency scales in the presence of finite damping at 
all frequencies. Numerical calculations based on this idea agree well with experimental data. 
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Introduction. Turbulence on the surface of fluids has 
been the focus of numerous experimental and theoret- 
ical investigations during the last few years [J, 0, IE 0, 
0, 0, H HI El ■ Interest in such phenomena, usually re- 
ferred as wave turbulence (WT), arises both from their 
great importance in terms of basic nonlinear physics 
and from numerous applications in engineering and the 
life sciences. One of the best known applications is in 
weather prediction, using information received from the 
measurements of the spectrum of waves on the sea sur- 
face. Turbulence in a system of capillary waves is of 
interest because its dynamics on this wavelength scale 
plays a significant role in the transfer and dissipation of 
energy on the liquid surface. Studies of nonstationary 
phenomena in systems of capillary waves are of partic- 
ular importance because they could provide direct in- 
formation about nonlinear wave interactions. In spite 
of the large number of experimental studies devoted to 
the nonlinear dynamics of surface waves, however, there 
are only a few recent reports 0, H, fl Hi of experimental 
studies of capillary turbulence that are directly compa- 
rable with theoretical predictions. 

Recently, the use of liquid hydrogen as a model 
medium has brought significant progress in the under- 
standing of capillary turbulence. The remarkable prop- 
erties of liquid hydrogen (its low viscosity and the high 
nonlinearity of its capillary waves, and the possibility 
of driving the electrically charged surface of the liq- 
uid directly) have allowed us to observe for the first 
time: the formation of the Kolmogorov power-law tur- 
bulent spectrum over a wide range of frequencies (10 2 - 
10 4 Hz) [4(; the modification of the scaling index of the 
turbulent spectrum in dependence on the spectral char- 
acteristics of the driving force [l(J; and the cut-off of 



the power spectrum of capillary turbulence at high fre- 
quencies due to the change of the energy transfer mech- 
anism from nonlinear waves transformation to viscous 
damping [ill ]. Many important properties of waves on 
the surface of liquid hydrogen are similar to those of 
conventional liquids like water (cf. the capillary lengths 
for liquid hydrogen A = 0.19 cm at T = 15 K, and for 
water A = 0.28 cm at T = 293 K), providing an addi- 
tional argument for using liquid hydrogen as a perfect 
test fluid for accurate tests of WT theory. The latter 
are of considerable importance, because WT theory is 
used to describe turbulent processes in a wide va riety 
of media, including e.g. plasmas [12], astrophysics [Ta j, 
ocean surfaces [13] , acoustic turbulence in superfluid He 
II [HI], and nonlinear optics [HI ]. 

In this Letter we report the first observations of decay 
of the turbulent state in a system of capillary waves on 
the surface of liquid hydrogen. We have observed that 
the decay of the stationary turbulent spectrum starts 
in the high frequency domain, with the energy remain- 
ing localised in the low frequency range of the turbulent 
spectrum, i.e. near the driving frequency. At low fre- 
quencies the turbulent spectrum remains close to its 
unperturbed shape for a relatively long period of time 
after the driving force is switched off. This observation 
differs significantly from what might be expected from 
the self-similar theory of nonstationary WT processes 
[l[, where the evolution of the spectrum is considered 
in the range of frequencies where viscous damping can 
be neglected. 

We show that where there is finite damping at all 
frequencies, viscous losses cause qualitative changes in 
the evolution of the turbulent spectrum after removal of 
the driving force: rather than a propagation of pertur- 
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bations from low to high frequencies, caused by the cas- 
cade transfer of energy (the scenario considered in [l| ) , 
there is a relatively fast decay of the high frequency do- 
main of the spectrum. Although our observations are 
new, and came as a surprise, we show that they can 
still be understood qualitatively within the framework 
of the general WT theory [l[ . 

I00 °1i * u\i i k i a (a) 

a. ItMMHHBMHHNNPFv ' r ' '' ■ ' 



1 1.5 ' 2.5 3 

t (sec) 




1 ' l!s ' 2 2.5 3 

t (sec) 



ID 1 ' 

10 

/I , 10 " 


. . ..i i ■ , , * 


2 T ^huu. 

< r , , 


100 1000 



100 



co/2ji (Hz) 

FIG. 1: (a) The measured signal P(t). The periodic driving 
force was switched off at time t = 1.8 s. (b) Evolution of 
the turbulent power spectrum during the decay, calculated 
over P(t). Grey shading indicates frequency components 
in the power spectrum whose P% exceeds the threshold 10 4 
(a.u.) in the graph below, (c) Instantaneous power spectra 
calculated at times indicated by the arrows in (b): curve 1 at 
the insert corresponds to time t = 2 s; curve 2 corresponds 
to i = 2.5 s. The spectra shown represent an ensemble 
average over 10 identical measurements. The dashed line in 
(c) corresponds to the power-like dependence P% ~ uj~ 7 ^ 2 
predicted by theory 0, 0] . 

Experimental observations. The experimental ar- 
rangements were similar to those used in our earlier 
studies of steady-state turbulence on the charged sur- 
face of liquid hydrogen [17| . The measurements were 
made using an optical cell inside a helium cryostat. Hy- 
drogen was condensed into a cup formed by a bottom 



capacitor plate and a guard ring 60 mm in diameter 
and 6 mm high. The layer of liquid was 6 mm thick. 
The top capacitor plate (a collector 60 mm in diame- 
ter) was located at a distance of 4 mm above the surface 
of the liquid. A two-dimensional positive charge layer 
was created just below the surface of the liquid with the 
aid of a radioactive plate placed at the bottom of the 
cup. The temperature of the liquid was 15-16 K. The 
waves on the charged surface were excited by a periodic 
driving voltage applied between the guard ring and the 
upper electrode. They were detected from the varia- 
tion of the total power P(t) of a laser beam reflected 
from the oscillating surface, which was measured with 
a photodetector, sampled with an analogue-to-digital 
converter, and stored in a computer. Given the size of 
the light spot, the correlation function I u = ( 1 77^ | 2 ) of 
the surface elevation r](r,t) in frequency representation 
is directly proportional to the squared modulus of the 
Fourier transform of the detected signal, I u = const P% 
at frequencies above 50 Hz. 

In systems of finite size, the power spectrum of cap- 
illary turbulence is discrete; but, at frequencies much 
higher than the lowest resonant frequency (which in 
our experiments was ~3 Hz), the intrinsic spectrum of 
resonant frequencies becomes quasi-continuous due to 
viscous broadening of the resonances and/or nonlinear 
broadening at pumping rates exceedinga critical value 
(see theory and observations [I II, i, S] ) . To es- 

tablish the steady turbulent state at the surface of the 
liquid, a 95 Hz ac driving voltage was applied for ~10 
s. It was then switched off, and we observed the re- 
laxation oscillations of the surface. The instantaneous 
power spectrum P% of nonstationary surface oscillations 
was calculated by using a short-time Fourier transform 
of the measured signal P(t) [lj|. 

It is clearly evident from Fig. 1 that, during decay of 
the turbulence, it is the high frequency components of 
the power spectrum that are damped first. The ampli- 
tude of the main peak at the driving frequency remains 
larger than the amplitudes of peaks at the harmonics, 
all the time, both before and after the driving force 
is switching off, i.e. the surface of the liquid continues 
to oscillate mainly at the driving frequency. The wave 
amplitude in the turbulent distribution decreases homo- 
geneously during the decay, and the power-law depen- 
dence of the spectrum persists, even for low frequencies 
down to 100 Hz. 

Numerical calculations. To try to understand the pe- 
culiar form of decay, we modelled numerically the time 
evolution of the initial steady-state turbulent power 
spectrum of the capillary waves. According to WT the- 
ory [J 0] the evolution of a turbulent spectrum can 
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be described by the integro-differential kinetic equation 
(KE) for the classical occupation numbers rik of sur- 
face waves, where k is a wave vector. In our calcula- 
tions we used the so-called "local approximation" for 
the KE, where the collision integral was approximated 
by a differential operator. This approach is similar to 
the well-known Fokker-Plank equation approach in the 
kinetic theory of gases [20], and is based on nonlinear 
interactions between waves of comparable frequencies. 
Such a differential analogue of the KE can be written in 
a unique way based on the known symmetry and scaling 
properties of the turbulent state. A similar phenomeno- 
logical model of WT has been used successfully before, 
e.g. in studies of gravitational waves on the surface of 
liquid [14] and for optical turbulence [3], yielding re- 
sults that are in a good agreement with the "exact" 
WT theory. To accord more accurately with the real 
conditions of the experiment, we introduce into (1) a 
viscous damping term for waves at all frequencies, in 
addition to the nonlinear term. Capillary turbulence 
has not hitherto been studied using this approach. 



Following the general recommendations of Refs. [14 . 
flij we write the corresponding equation for the envelope 
curve of the turbulent spectrum as: 
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For convenience in comparing numerical results with 
our experimental data, we use the occupation num- 
ber in w-representation, n u (t) — nk(t)\k = k(uj) a s the 
main characteristic of the turbulent distribution, where 
k — k{oj) is the function inverse to the dispersion law for 
the capillary waves u> = lj^ ■ The first term on the right 
hand side of ([T]) plays the role of the collision integral 
in the KE, 7^ oc oj 4 / 3 is the damping coefficient of cap- 
illary waves, the term f u (t) models the direct action of 
the driving force, and C is a dimensional constant. For 
capillary waves the correlation function I u = const n w 
has the same scaling properties In our simula- 

tions we renormalised time by a dimensional constant, 
t — > t/C, and renormalised frequency, w — > Cui. We 
have assumed that the renormalised driving frequency 
ujd/2-K = 100 to facilitate comparison of our numerical 
modelling results for the free decay shown in Fig. 2 with 
the experimental data of Fig. 1. The full curve shows 
the envelope of the initial steady-state turbulent distri- 
bution n u . For the frequency range 200 < oj/2tt < 800 
it can evidently be approximated by a power-like de- 
pendence riuj ~ w -7 / 2 , in good agreement both with 
the existing theory of capillary turbulence 0, 0] and 
with our current and previous experimental observa- 
tions [l(| • The cut-off in the power spectrum at the fre- 
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FIG. 2: Decay of the turbulent spectrum. The occupation 
number a function of renormalised frequency calcu- 

lated using (fl| at three renormalised times: t = 0, the initial 
steady-state spectrum (full curve); t — 3 (dashes); t — 10 
(dots). The straight line corresponds to ~ uj~ 7 ^ 2 as pre- 
dicted by the theory of capillary turbulence [1, 0] ■ 



quency uj^jl-K m 800 is caused by viscosity. The dashed 
and dotted curves show the evolution of the envelope 
after the driving force was removed. The destruction of 
the power-like dependence clearly begins from the high 
frequency domain of the turbulent spectrum. The en- 
ergy of the oscillations is concentrated mainly in the low 
frequency spectral domain. That there is no destruc- 
tive front propagating from low to high frequencies is 
entirely consistent with our experimental results. 

These results may be explained qualitatively in terms 
of the general theory of WT if we assume a fast redis- 
tribution of energy between different frequency scales 
inside the inertial range, leading to suppression of re- 
laxation processes. Such a redistribution will stabilize 
the power-like dependence of the turbulent spectrum 
at low frequencies. A similar evolution of the turbu- 
lent spectrum was observed in [ill ] with a smooth de- 
crease of driving amplitude, where the turbulent sys- 
tem remained continuously in its steady state. Based 
on these observations we can claim that the evolution 
of the freely decaying turbulent spectrum has a quasi- 
adiabatic character. The kinetic time of nonlinear wave 
interactions Tk(u>) plays the role of a fast time; the time 
t v (u)), characterizing viscous damping of the waves, 
plays the role of a slow time within the inertial fre- 
quency range. From ([T]) one can estimate the frequency 
dependences of Tk(uj) oc oj~ 7 / 6 and t v (u>) oc w~ 4 / 3 , and 
of their ratio as r(u>) = t^{oS) / t v {oj) oc (aj/wt,) 1 / 6 . It 
was shown in [11] that the inertial range is limited at 
high frequencies by the condition r(u>b) ~ 1, and that 
inside the inertial range the parameter r is small, i.e. 
nonlinear wave interactions occur faster than the vis- 
cous damping processes. At these frequencies the ratio 



4 



r(u>) can be considered as a small adiabatic parameter. 
At frequencies of the order of Wb, r(uJb) ~ 1, so that 
the nonlinear and viscous processes have comparable 
rates, and the adiabatic condition is violated. During 
the decay process the amplitudes of the waves decrease, 
so the kinetic time decreases as well 111! . Hence 



the shift of the boundary frequency u>b towards the low 
frequency domain. At large times after the beginning 
of the decay, u>f, becomes comparable with the driving 
frequency cu d . Dissipation is then starting to play an 
important role even at low frequencies. Thus we can 
claim that, for the understanding of nonstationary tur- 
bulent processes, viscous losses in a turbulent system 
are of central importance at all frequencies above the 
driving frequency: we have seen, both in experiments 
and in computations, that finite damping of the waves 
changes qualitatively the character of the turbulent de- 
cay, compared to the scenario envisaged earlier. Note 
that these considerations do not apply to stationary tur- 
bulent phenomena, where dissipation can be neglected 
completely over a wide range of frequencies 0, (9| . 

Conclusion. Our experimental observations and nu- 
merical calculations have shown that the decay of capil- 
lary turbulence on the surface of a normal liquid begins 
with the damping of high frequency waves, so that the 
high frequency part of the power-like turbulent spec- 
trum is destroyed first. The energy-containing range of 
frequencies does not shift toward the high frequencies 
in practice, but remains at frequencies of the order of 
the driving frequency. This remains true at all times 
during the decay, even at large times when the turbu- 
lent regime is replaced by the viscous damping of the 
waves. Transition processes in the turbulent system are 
damped by the fast nonlinear redistribution of energy 
between waves whose frequencies are inside the iner- 
tial range. These general conclusions may be useful for 
future studies of nonstationary turbulent processes in 
a wide class of problems, e.g. in hydrodynamics, non- 
linear sound in solids and liquids, astrophysics, and 
plasma physics. In particular, our results support the 
phenomenological propositions made in [2l| to account 
for the observed long-time properties of vorticity in the 
decay of quantum turbulence in superfluid helium. It 
could be also useful in relation to recent studies [l8| 
of damping of the monochromatic capillary wave on a 
liquid surface. 
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